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FOREWORD 


The program reported herein was performed by the General 
Electric-Space Division, Valley Forge, Pa., for the George C. Marshall 
Space Flight Center, Marshall Space Flight Center, Alabama, under 
Contract NAS 8-28167. The performance period for the work was 
2 December 1971 to 29 June 1973. The principal investigator was 
Edward J. Kuhar and the program manager was Clyde V. Stahle. The 
NASA Technical Monitor was Dr. John R. Admire who provided valuable 
guidance throughout the course of the program. 

The results of the study are described in the main volume of this 
report and include the theoretical development of the dynamic transforma- 
tion method, numerical results from the application of the method to 
several sample problems, and some comparisons with other available 
methods of analysis. 

The separate addendum to this report provides the user instructions 
for the DAMUS computer program (Dynamic-transformation Adapted to 
Modal-synthesis Using Stiffness Coupling) v/hich implements the method 
developed under this program. 
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Nomenclature 



mass matrix for substructure in ^x^ physical 
coordinates 


M 


stiffness matrix for substructure in £xj physical 
coordinates 



matrix of substructure eigenvectors in £x 
coordinates 


K] 

[ k t ] 


mass matrix for total structure in jx| physical 
coordinate s 

stiffness matrix for total structure in ^xj physical 
coordinates 

matrix of eigenvectors for total structure in ^x^ 
physical coordinates 



generalized mass matrix for total structure in 
modal coordinates 



generalized stiffness matrix for total structure in £qj 
modal coordinates 



matrix of eigenvectors for total structure in £qj 
modal coordinates 


[t] 

H 

H 


dynamic transformation matrix 

transformation matrix defining the relationship of the 
reduced coordinates, £q^j » to the kept coordinates ^q^ 

physical coordinates 

generalized modal coordinates 

system circular frequency 


U) 


substructure circular frequency 


P 


reduction circular frequency 


n 


number of degrees of freedom in structure 
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Subscripts 


A 

CPL 

k 

r 


refers to the i subset, term, or substructure 
incremental mass from coupling spring 
incremental stiffness from coupling spring 
kept coordinates 
reduced coordinates 


k 

r 

A 

I 


ij 


Superscripts 

kept coordinates 
reduced coordinates 
attachment coordinates 

interior coordinates not attached to any other sub 
structure 

revised value 

second time derivative 

particular submatrix partition 
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SECTION 1 


INTRODUCTION 

This document describes a Fortran IV computer program used to per- 
form modal synthesis by stiffness coupling using the dynamic transformation 
method. The program has been named DAMUS (Dynamic-transformation 
Adapted to Modal- synthe sis Using Stiffness-coupling). The program begins with 
the entry of a substructure's mass and stiffness matrix. The eigenproblem for 
the individual substructure is solved. Provisions are included for a maximum 
of 20 substructures (100 DOF max/substructure) which may be coupled by 100 
stiffness matrix springs (100 DOF/spring). The substructures are then coupled 
together via coupling springs, and the dynamic transformation is used to reduce 
the size of the eigenproblem. After solving for the coupled system eigenvalues 
and vectors, the user may elect to backsubstitute selected modes. The total 
number of modes treated by the program is 300 consisting of 100 kept coordin- 
ates (maximum eigenvalue size) and 200 coordinates reduced by the dynamic 
transformation. For user flexibility, six major entry points have been included 
in DAMUS. 

Input data for DAMUS is mainly accomplished by the READ and READIM 
FORMA subroutiites. Output data to be saved is written on files generated by the 
WTAPDS and WTAPSS subroutines written specifically for DAMUS. Those files 
which contain data to be saved should be copied to tape after the execution of 
DAMUS. A total of 12 files have been defined for use by the program. Depending 
on user options, the number of files used at any one time will vary; and at no 
time will all twelve be used simultaneously. 
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SECTION 2 

THEORETICAL DISCUSSION 


2. 1 BASIC THEORY FOR STIFFNESS COUPLING 

The stiffness coupling method of modal synthesis assembles the com- 
plete structure in the same manner as the displacement method for structural 
analysis. The total structure may be represented by a number of substructures 
connected through flexible links. Each substructure is analyzed without the 
flexible links to determine the component vibration modes with free attach- 
ment coordinates. The flexible links are represented by a stiffness matrix 
relating the interface forces from one set of substructure attachment coordin- 
ates to another. 

The method of substructuring for stiffness coupling may best be illustrated 
by considering a total structure consisting of only two substructures. The 
general undamped equation of motion for the ith substructure in terms of its 
generalized mass matrix, QvW^D, and generalized stiffness matrix, aa , 
is given by 


4- (1) 

where the coordinates describe physical motions of the mass points. 

Each substructure has two sets of coordinates which will be referred to as 
attachment coordinates, £ "X-i j , and internal coordinates, { . The 

attachment coordinates are those degrees of freedom (DOF) which are connected 
to another substructure via a stiffness matrix or coupling spring. The internal 
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coordinates are those DOF which are not connected to any other substructure. 
It is important to note that for stiffness coupling the coordinates belonging to 
one particular substructure are not common to any other. If Y\ - L represents 
the size of the ith substructure and Vv the size of the total structure com- 
prised of substructures, then V\ will be given by 


V\ = 


A 




( 2 ) 


Having defined two sets of coordinates for each substructure, Eq. (1) may be 
written in partitioned form as 


It I IX 

i 

21 I *X 

I Wi 

• I _ 

y\i* y\\. 



o 


( 3 ) 


Now consider the total structure to be described by a mass matrix, [ > 

and a stiffness matrix, LKtI , such that 

[Mt]{x t ] + [Kt]{*t] =■ o (4) 

nnn rv*| **1 


where 


[Mr] 

v^. x \n 


• I I 

M t j^O 



*i 



x I 


( 5 ) 
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If we describe the connecting structure between £ 7C* | and 

by a mass matrix, , where 

I I » 

*. A 

'' [ 7 , . A 

(h.Vn.XW'.^jL ! Al J 1 


n, 


and a stiffness matrix, C*c pt ]. relating the nodal forces, £ ^ , to 


by the equation 


c;c u,>' 






(n, *♦*,*> I 


u A . A 
n*. 


Then [Mr] and W may be written in partitioned form using the submatrices 


defined by Eqs, (3), (6) and (7): 


M ■ 


l' " I It I 
*V\, + «na, | W\, | 0 


I l II It 


zi 


Z-3 
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and recognizing the partitions belonging to each of the substructures, Eq. (4) 
may be written as 


Yr\, \^0 

O I 


Ha 

V\*h 


'CPL 


W 


v 

— -> - o 


( 10 ) 


V\*h *»* 


•«> *\ Cr i 

i 

’ + ' ~ ~ 1 "7~ 

^2, 0 | 

J 

V“\*l V\*V-\ 

For each substructure defined by Eq. (1), a set of eigenvalues, OV) , 
and a set of mass normalized eigenvectors, 03 , can be obtained such that 

~ T Tw,,i r - r ti 

(Ha) 


03 MM = [i] 
[<i»0 T [-«0t*i] ^ Ovj 


(lib) 


Using the results from Eq. (11), we can now express Eq. (10) in terms of a 

The coordinate transformation 


set of generalized modal coordinates, W 
is given by 

{*} = ! 

v\* v\ n* I 


where 


{%} «[♦]{»} 

Y\* I V\ 

W • fel 


r»«u 


w-(£l 


» 


w ■ [-Hi] 

Y\A V> 


( 12 ) 

(13a) 

(13b) 


Substituting the transformation for into Eq. (10) and premultiplying by 

[« yields 


{[i>W T [*vlW}{fl 


F* • ^ 1 

20, I 0 

l il 




0 (14) 


h»h 
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Solution of Eq. (14) will result in a set of eigenvalues, , and then 

corresponding eigenvectors, DO Because of the similarity transforma- 
tion used to obtain Eq. (14) from Eq. (10), the eigenvalues of Eq. (10) will 
be equal to those of Eq. (14) and the corresponding eigenvectors, , 

of Eq. (10) will be given by 

ha = mi>] 


(15) 






Eq. (14) represents the most general form of the equation of motion for stiffness 
coupling. This equation is generally solved by partitioning the M coordinates 
into two groups, kept and truncated. The truncated coordinates correspond to 
the high frequency substructure modes and are completely omitted from the 
equation of motion. Those degrees of freedom remaining, the partitioned set 
of kept coordinates, determine the final reduced size of the eigenvalue problem 
to be solved. 


The general form of Eq. (14) may be simplified further by including the 

r 

correct partitions from Eq. (6) in each corresponding at the 

substructure level. This is reasonable because we are assuming that there is 
no inertial coupling between substructures. This will result in E = 0 and 

Eq. (14) reduces to 

CO{fr} ♦£*]{»}-<> (i6) 

V\*V\ WKl V\«« V\*t 


where 


DO- 


R x i - 

o 

-o-V 


r>*y\ 


w’ww 


(17) 
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After solving for the Cu^'s and 4>{ 's, the only lengthy calculation left to 


be performed in order to obtain Eq. (16) is the matrix triple-product involving 
F Kepi. J . If we partition each row-wise in terms of its and 

coordinates such that 


hi xh;. 

Eq. (13b) can be written as 




O] = <p* 


K *,* 

‘P , 1 


By forming the triple product using Eqs. (9b) and (19) and 

refactoring the result in terms of the matrix partitions, the resultant form of 


the triple product may be expressed as 

I> A ] T [^] [+*] 


(20a) 


■ 

cn = 

. ° ! < 

(h, A 4h 




M ~ 

"a » ' 
^CPl ^ 

^1' 1 
Kcpl j - 

u. A 

(n,*+h l V(«,''+h*) 


(20b) 


(20c) 


and the final form of M in Eq. ( 17) may be expressed as 

DO + Cn T [4.][^] 


h* n 


o iw a 


W\x v\ 


2-7 



2.2 DYNAMIC TRANSFORMATION 


As a result of omitting the higher substructure modes, the solutions 

from the truncated Eq. (16) will have errors introduced. The truncation errors 

can be greatly diminished by including the modes that would have been truncated 

through a dynamic transformation. Instead of truncating or omitting modes, 

all modes can be included through a transformation that relates the "reduced" 

modes not contained explicitly in the solution to the modes that are "kept. " If 

2 

si i corresponds to an exact eigenvalue of Eq. (16), the relationship between 
the eigenvalue and its eigenvector may be expressed in terms of the kept, ij-i' 

, coordinates as: 


and reduced, 


W'" 

k^j If 


( 22 ) 


V\*W l 


V\ * Y\ 


Y\* t 


where the m corresponds to those modes previously truncated. If we 
designate Y \ £ as the total number of modes kept from all the substructures 
and Vhy. as the total number of modes reduced, then 


Y\ = V)fc +• VXy 


(23) 


2 ^ 4- 

Expanding Eq. (22) into two equations for some general frequency, p = , 


yields 




(24a) 


I 






t'W - * mb' 


(24b) 


i 


D**! tty* tty VNyX) 
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Solving Eq. (24b) for If! in terms of give! 

V\y*l ^*1 


( 25 ) 


where 


W i } 1 L< y *] 

V\yX.V\y V\yXY\ h 


( 26 ) 


Using Eq. (25) for some ’’reduction frequency", p , we can write 


(a* Co* 

\±.i = .1 v = 

if\ h 

h * I h * I 


I 

~R 


(27) 


V\**l 

The dynamic transformation matrix, Lt] , is then defined as 


r 1 1 

m - b'J 


(28) 


The reduced equation of motion is obtained directly by substituting the coord- 
inate transformation 




(29) 


into Eq. (16) and pre-multiplying by the transpose of [TJ . The reduced 
generalized mass and stiffness matrices can be written in the partitioned forms 
given by Eqs. (22) and (28): 

CM*] = [I] + [R] r [R] 

by-*** 


(30a) 


D<*] = [<“] + 2 + [R] T [K-J[R] 

" t ‘" 6 "»* h ' **'•'* ^ 


(30b) 
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Conventional methods of determining eigenvalues may be applied to the 
reduced equation of motion to obtain a set of eigenvalues, , and a 

corresponding set of mass normalized eigenvectors, [y*J From the coord- 
inate relationship defined by Eq. (25), the reduced eigenvectors, O r ] , 
corresponding to the hi reduced coordinates are given by 

O v l = M tV] 


where 




r 

L ^ J ^ tty 


and the physical eigenvectors for the total solution will be given by 




h*vi£ 

■A 

This solution will be exact for any which is the same as the reduction 

frequency, p , used in developing [T] • 

Significant improvement can be obtained in the modes and frequencies by 
applying the Rayleigh-Ritz method in conjunction with the dynamic transforma- 
tion for individual modes. This part of the dynamic transformation will be 

A 

referred to as backsubstitution. For each to be considered, a revised 

mode shape, [ ¥<; J, can be determined by substituting in C8] 

and mass normalizing the mode shape: 
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(34a) 


i>n = i>i][vi*] 

- t^] T rC] 


(34b) 


NFi. = '/VoT 


(34c) 


lv?] = nf £ [ x«] 


(34d) 


where N is the normalization factor used to mass normalize one for 

£ 

some SLi 


A new estimate of 



will be given by 


Another measure of the accuracy of the solution is provided by the change in 
the eigenvalues and is provided by the eigenvalue ratio defined as 


& /-o 


Z 
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SECTION 3 


PROGRAM ENTRY POINTS 

In order to provide user flexibility and at the same time minimize 
computation of basic data changes, six major entry points have been 
established for DAMUS: 

1. Basic substructure data entered, subsystem eigensolutions. 

2. Coupling spring stiffness data entered and stiffness contribu- 
tions calculated. 

3. Selection of modes to be kept/reduced/truncated, generalized 
mass and stiffness matrices calculated. 

2 

4. p value for dynamic transformation entered, eigensolution 
for system using the dynamic transformation. 

5. Calculation of physical eigenvectors for total system. 

6. Backsubstitution of selected modes from system solution. 

Entry into the program is accomplished by designating a specific entry 
point (EP). Termination of the program is accomplished by designating 
the last EP the user desires to execute. 

The Fig- 3-1 flow chart is included to show the general flow of 
the program. A corresponding flow chart showing the main Fortran 
subroutine called by the program may be found in Appendix B. The flow 
chart for input/output files required by the program at each EP is in 
Section 4. This flow chart corresponds in form to the one shown in 
Fig. 3-1. 
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Read p transformation 
frequency; reduce mass and 
stiffness matrices; obtain 























A special vector input notation is used for reading in vectors used 


to select DOF orders. When it is required for the user to select a group 

of DOF's to be re-ordered or to designate a group of coordinates to be 

printed from a substructure, a vector of identifying DOF's is read into 

DAMUS by the READIM FORMA subroutine. The order in which the DOF 

numbers appear determines the particular sequence of DOF desired. If 

a sequential set of numbers from N1 to N2 is desired, the input may be 

abbreviated by the user. Inclusive groups of numbers to be generated in 

ascending order may be specified at any one time by using three elements 

of the input vector IV where 

IV(I) = N 1 
IV(I+1) = 0 
IV(I+2) = N2. 

The integers from N1 to N2 will be sequentially expanded in the IV matrix 
starting from the IV(I) location. The last element of IV must be negative 
if this abbreviated form of input is used. 

For an example, let us consider a substructure with 20 DOF's. It is 
desired to print only 10 of the substructure DOF's in a different order. The 
vector for reordering the DOF's would be designated as a 1 x 10 on the 
READIM header card. If the order to be printed is given as 
7, 8, 9, 10, 2, 1, 15, 16, 17, 18 
then the shortened input vector would be given by 
7, 0, 10, 2, 1, 15, 0, - 18. 
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If all the DOF's were desired to be printed in their original order, the 
input vector (1X20) would be 
1 , 0 , -20 

and the expanded vector of numbers 1 through 20 would be generated by 
the program. 


Entry Point 1 

The mass and stiffness matrices, [mil, C4 t l , for each substruc- 
ture are read into DAMUS at EP-1. The substructures are defined in the 
program by a user supplied number which ranges from 1 to 20. Since 
each substructure is to be identified in the program by a distinct number, 
their input may be in any order. Input data for the mass matrix may be in 

two forms: an n. x n square matrix or a 1 x n. row vector. Before 
i i i 

solving the substructure eigenproblem, the attachment coordinates must 

A 

be identified and partitioned into the x. set. An input vector IDDOF is 

A 

used here to specify the N,, coordinates. Since the special input notation 

as previously described is to be used, only the n. coordinates need be 

specified. The program will complete the vector for re-arranging the mass 

A 

and stiffness matrices. The n. coordinates must be partitioned in the 

same order as the coupling- spring DOF's. If one substructure couples 

to several others, then each set of attachment DOF's must be specified 

in the order in which they will be used. For example, assume that sub- 

A 


structure i couples to 3 other substructures. The "X. coordinates may 
then be partitioned as 
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Al 


AZ 


Al 


hi 


vv 


AZ 




Xi 


Al 


K 


.Al 


Al A2 

(If X. = X. , then only 2 partitions need to be specified). For EP-2, 
it will be necessary to input the n. * locations since each spring only 
couples 2 substructures at a time. Thus, there will be at most, 3 sets 

, a 

of Cp. associated with substructure i. After re-arranging the sub- 
structure DOF's, the complete vector will be printed out to allow the user 
to identify the substructure mode shapes. An additional input vector is 
required here to select which substructure DOF's are to be printed after 

the system is coupled. The DOF order desired refers to the original rm, 

A i 

k., read into core, not the partitioned X. , X. coordinate set. If one 

desired all of the coordinates to be printed, the abbreviated vector 

( 1 x n.) would be given by 1, 0, -m in the READIM format. An option is 

included at this point to rotate the m., k into system coordinates by 

i i 

reading in a 3 x 3 direction cosine matrix. The substructure may not 
be rotated if its size is not divisible by three. 

Another option for EP-1 is one for altering or adding substructure 
data. The user may add new substructures at any time. If a substructure 
or group of substructures need to be altered, only those to be changed can 
be entered. The program tapes will be updated to reflect the changes 
without regenerating previous unaltered substructural data. 
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Entry Point 2 


Coupling spring stiffness data is entered here. Each coupling 
spring stiffness matrix is associated with only 2 substructures. The 
user must specify which two substructures are being coupled, the number 
of attachment DOF's associated with each substructure, and the starting 
location of the n. A 'y 

coordinates that were specified by the IDDOF vector 
in EP-1. The program then performs the product 

Dt> A j T r 

and saves the result for assembly in EP-3 where the ^ contains only 
Ai 

those %- coordinates corresponding to the Kcpl;_ DOF's. As in EP-1, 
the user may add additional coupling springs corresponding to more sub- 
structures added. If a substructure was altered, this EP must be executed 
again to reflect the substructures' altered eigenvector. Only those coupling 
springs of direct concern need be calculated. The program will update 
the tapes without regenerating previous unaltered substructure data. 
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Entry Point 3 


At this EP, the user specifies which modes are to be kept and re- 
duced. The input matrix KEEP (2 x MAXSQB) is used where the column 
corresponds to a substructure. The ( 1 , i) location specifies the total 
number of lowest modes kept for substructure i and the (2, i) row specifies 
the total number of next highest modes to be reduced. All other modes 
will be truncated. The restrictions for the total kept and reduced modes 
are given by 

MAxsua 

2 KEEP (l» ^ too 

£=l 

MM.SU0 

K£EP (Zji) 200 

<-■* I 

The generalized stiffness matrix is then assembled in its kept and reduced 
partitions. There must be at least 2 modes kept from each substructure 
read into DAMUS. For KEEP(1, i) = 10 and KEEP (2, i) = 20, the program 
will select modes 1 to 10 for the kept partition from substructure i and 
will place modes 1 1 to 30 in the reduced partition. 

Entry Point 4 

The reduction frequency, p , is entered at EP-4, and the dynamic 
transformation is applied to obtain the reduced mass and stiffness matrices. 
The eigenvalues for the coupled substructures are now obtained. 
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Entry Point 5 


This entry point calculates the coupled physical eigenvectors* 
which are printed in substructure groups. The order in which they are 
to be printed is determined by the order the substructures were originally 
read into DAMUS. The user must specify how many modes are to be 
printed by specifying the first and last mode number of the group of 
modes desired. Considerable computer time may be saved by selecting 
coordinates and printing a few of the modes since only those selected will 
be calculated. The user may desire to re-enter the program to select 
new coordinates to be printed or more modes. If the option to obtain new 
coordinates is selected, the user must specify which substructures will 
be printed and a new identification vector for selecting the coordinates 
for each substructure being considered. When selecting only different 
modes to be printed, those coordinates previously defined in EP-1 will 
be used to calculate physicals. An option also exists for making the co- 
ordinate selection a permanent change and will cause the basic program 
tapes to be updated to reflect any new ordering specified by this entry 
point. The order in which the substructure groups will be printed may 
be altered this way. 


^Referred to as "physicals" for convenience. 
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Entry Point 6 


This EP is for backsubstitution. The user may select in any order 
up to 30 modes at a time to be backsubstituted. The same physicals as 
specified in EP-5 will be calculated and printed. If the program is entered 
here, it will use the system tapes previously generated and saved to cal- 
culate the physicals. For new coordinate selection in EP-5, the correct 
save tape must be used to obtain the physicals. If EP-6 is executed after 
EP-5 and a new coordinate selection was made, the physicals calculated 
will be those currently specified in EP-5. The new eigenvalues will be 

printed along with the normalization factor used to mass normalize the 

— A / fe. 2 

vectors and the frequency ratio (-TL^ / S\ -i ) used as a measure of the 
frequency change. 
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SECTION 4 


INPUT/OUTPUT DATA 


4. 1 INPUT DATA 

This section describes the necessary input data required to 
execute DAMUS and is ordered by entry points. The data cards and 
the variables appearing on them are listed in sequential order with 
the input format or FORMA subroutine specified. The definitions for 
the variables specified on the cards are given to clarify their meaning. 
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Input Data for Starting Program 


Card 

No, 

1 

2 

3 

4 

5 

6 


IRUNNO 
UNAME 
TITLE I 
TITLE 2 
FOD 

NTAPE(I) 


IENTR 


IEXIT 


Input Order 


F ormat 


IRUNNO, UNAME 
TITLE 1 
TITLE 2 
FOD 

NT A PE ( 1 ) , NTAPE(2) 

IENTR, IEXIT 


(A6, 4X, 3A6) 
(12A6) 

(12A6) 

(E10. 0) 

NT A PE (12) (1215) 

(215) 


Definitions 


Run No. 

User's Name 
Fir st T itle 
Second Title 

Final off-diagonal value for diagonalizing a 
matrix M using the method of Jacobi. 

Twelve tape and/or file units used by the program 
and assigned by the user. Only those tapes 
actually used must be specified. 

Entry Point for entering program. There are 6 
entry points in all defined by the integers from 
1 to 6. If the program is being entered at Entry 
Point 1 other than the first time, then IENTR = - 1 
must be specified. 

Exit Point for terminating the program. The exit 
points are identical to the entry points with program 
termination occurring after the execution of the entry 
point specified by IEXIT, 
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4. 1. 1 Input Data - Entry Point 1 


Card 

No. Input Order Format 


1 

2 

3 

4 

5 



NSUBS, IOP (215) 

ISUB, KEPMOD, MOPT, IROT (415) 
IDDOF (1XN1) READIM 

K (NXN) READ 

M (NXN) or ( 1XN) READ 

R COS (3X3) READ 

IPDOF ( IXN2) READIM 


Cards 2 through 7 are repeated for each substructure defined. 


Definitions 


NSUBS 

IOP 


ISUB 

KEPMOD 


Number of substructures to be read into program 
at this time. There will be NSUBS sets of cards 
from 2-7 following Card 1. 

0 denotes the first time the program is entered 
at Entry Point 1 

1 denotes more substructures are to be added 
to a set of substructures previously defined 

2 denotes one or more substructures from a 
previously defined set are to be changed. 

Substructure identification number which ranges 
from 1 to 20. Each substructure must be identified 
by a different no. 

Total number of modes to be saved on tape for sub- 
structure "ISUB. " This number must not be less 
than the sum of the number of modes to "kept" plus 
the number of modes to be "reduced" for this parti- 
cular substructure. 


*( ) denotes optional input data and is not required if option is not exercised. 



MOPT 

TROT 


IDDOF 


K 

M 

RCOS 


IPDOF 


1 Mass matrix is square NXN 

2 Diagonal mass matrix to be read in as a 
1XN vector. 

0 Substructure already in system coordinates. 
Delete Card 6 from input. 

1 Rotate substructure by direction cosines, 
a 3X3 matrix, read in on Card 6. 

Special 1XN1 integer vector used to place connecting 
coordinates in first N1 DOF locations for "ISUB. " 
Order of connecting DOF's must correspond to 
those used in coupling spring. The order in which 
the DOF's appear in IDDOF determines the re- 
arranged DOF order. N1 = the number of elements 
in IDDOF. 

Input Stiffness Matrix, for "ISUB". (NXN) 

Input Mass Matrix for "ISUB". (NXN) or (1XN) 
Optional Direction Cosines Matrix (3X3). Total 
number of substructure DOF's must be divisible 
by 3 in order to exercise this option. 

Special 1XN2 integer vector used to select those 
physical coordinates to be calculated from "ISUB". 
The order in which the DOF's appear in IPDOF will 
determine the order in which they will be printed 
for substructure "ISUB." If no physical coordinates 
from "ISUB" are to be calculated, IPDOF (1) must 
equal zero. N2 equals the number of elements in 
IPDOF. All substructure DOF numbers refer to 
the original order in which they were read into the 
computer. 


4-4 



4. 1. 2 Input Data for Entry Point 2 


Card 

No. 

1 

2 

3 

NCPLS 

IOP 

ISUB1 

NROW1 

NS1 

ISUB2 

NROW2 

NS2 

KCPL 


Input Order Format 


NCPLS, IOP (215) 

ISUB1, NROW1, NS1, ISUB2, (615) 

NROW2, NS 2 

KCPL (NXN) READ 


Definitions 


Number of coupling springs to be read into 
program at this time. There will be NCPLS 
sets of cards 2-3 following card 1. 

0 denotes the first time the program is entered 
at Entry Point 2. 

1 denotes more coupling springs are to be added 
to a set of data previously generated 

2 denotes one or more coupling springs from a 
previously defined set are to be permanently 
changed. 

First substructure coupled by spring (ISUB no. ) 
Number of connecting DOF's in ISUB1 for this 
particular spring. 

Starting DOF location of coordinates in ISUB1 which 
were ordered by IDDOF. The substructure may 
have more than one set of connecting coordinates. 
NS1 defines the reordered starting location for a 
connecting set of DOF's. 

Second substructure coupled by spring 
ISUB2 must always be greater than ISUB1. 

Number of connecting DOF's in ISUB2 for this 
particular spring. 

Starting DOF location of coordinates in ISUB2 which 
were ordered by IDDOF. 

Coupling spring stiffness matrix for coupling ISUBl 
to ISUB2. DOF's for coupling spring are in same 
order as specified for the ISUBl, ISUB2 coordinates 
by IDDOF with the ISUBl coordinates appearing 
first. N is the size of the total number of coord- 
inates U9ed from each substructure. 
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4. 1. 3 Input Data for Entry Point 3 


Card 

No. Input Order Format 

1 KEEP (2XMXSUB) READIM 


Definitions 

KEEP ' = INTEGER Matrix defining how many modes are 

to be kept and reduced for all of the substructures. 
The KEEP(1,I) row defines the number of low 
modes to be kept for substructure I. The KEEP 
(2,1) row defines the number of modes to be reduced 
for substructure I. Modes from the substructures 
are arranged in ascending order by frequency. For 
some substructure ISUB, the first KEEP(1,ISUB) 
modes will be kept and the next KEEP(2, ISUB) modes 
will be reduced by the dynamic transformation. 

MXSUB is the highest numbered substructure read 
into the program. A zero in some KEEP(2, ISUB) 
location will include all modes in the kept set. 
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4. 1.4 

Input Data for Entry Point 4 


Card 

No. 

Input Order 

F or mat 

1 

LAMDAO (1X1) 

READ 


Definitions 


LAMDAO 


value used for the dynamic transformation 
for reduction. 
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4. 1. 5 Input Data for Entry Point 5 


Card 

No. Input Data Format 


I 

( 2 )* 

(3)* 


MD1, MD2, IOP 
ID ( 1XN1) 
IPDOF ( 1XN2) 


(315) 

READIM 

READIM 


There will be N1 Card (3)'s. 


Definitions 


MD1 

MD2 

IOP 


ID 

IPDOF 


Mode number of first mode to be printed. 

Mode number of last mode to be printed. 

0 no optional cards needed 

1 requires cards (2) and (3), but data files 
not updated 

2 requires cards (2), (3) with permanent update 
of files. 

An integer vector containing N1 substructure ID's 
for which new physical DOF's are to be defined. 

Ordering vector for selecting DOF's from each 
substructure defined by ID. There will be N1 
of these cards. 


*( ) Indicates optional data and is not required if option not exercised. 
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4. 1. 6 

Card 

No. 

1 

2 


NBKSB 


Input Data for Entry Point 6 


Input Data 


F ormat 


NBKSB 15 

IMODE ( 1XNMD) READIM 


Definitions 


Number of groups of modes to be backsubstituted. 
30 modes may be specified to be backsubstituted 
at any one time. NBKSB tells the program how 
many groups of 30 modes are to be selected. 


IMODE 


Integer vector selecting modes to be backsubstituted 
with NMD £ 30. There will be NBKSB card (2)'s. 



4.2 OUTPUT FILES/TAPES 

This section describes the data saved on each logical file. Figure 
4-1 presents a Tape Flow Chart indicating which files will be required 
for executing the desired entry points. The program is written to handle 
files instead of tapes. A computer utility program should be used to save 
the necessary files on tape for re-entry into DAMUS. All eigenvalues and 
eigenvectors are output by subroutine WTAPDS and may be read by 
RTAPDS for use in another program. 
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Figure 4-1. Tape Flow Chart 
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4.2.1 Description of Output Files/Tapes (Ref. Fig. 4-1) 


NTAPE(I) 

I 


Description 


1 


2 


3 


4 


5 


This tape contains basic data necessary for program 
logic. It includes substructure sizes, ordering data, 
coupling indices, and data locations on the various tapes. 
It is required for every Entry Point (EP) and will be up- 
dated to reflect any changes due to option selections. It 
consists of one logical record written in binary by the 
standard Fortran WRITE routine. 

This tape contains the eigenvalues and eigenvectors for 
each substructure. In addition, it has the identification 
vectors used to re-order the substructure DOF's for the 
initial solution and for printing physicals. This tape 
is output data from EP-1 and is used by EP-2, 3. If 
IOP = 2 for EP-5, this tape is needed as input. A new 
NTAPE(2) will be output on NTAPE(12) for subsequent 
runs. 

This tape contains the partitioned set of substructure 
eigenvectors to be used in calculating physical vectors. 
Those DOF's to be calculated are packed into 100 DOF 
blocks and saved in the order read in. This is output 
from EP-1 and is required as input for EP-5 and EP-6, 
only if IOP = 0 for EP-5. 

This tape contains the partitioned coupling spring data 
used for assembling the generalized stiffness matrix. 

The substructures are coupled 2 at a time. This tape is 
output from EP-2 and is required input for EP-3. 

This tape contains the assembled generalized stiffness 
matrix in partitioned form plus the eigenvalues and 
participation factors from the system solution. Data is 
output from EP-3 and EP-4 on this tape and is required 
as input for EP-5 and EP-6. Re-entry into the program 
at EP's 3 and 4 will destroy previously generated data 
for that particular EP. If it is desired to save the pre- 
vious data for some future reference, a new tape should 
be used; i. e. , copy old tape data to new tape and use if 
data from EP-3 is needed. 
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NTAPE(I) 


Description 


I 

6 


7 

8 


9 


10 


11 


12 


This is a scratch file. It is basically required only 
for EP-4. For particular options, it may be required 
for EP's 1, 2, 5, and 6. If IOP > 0 in EP-5, NTAPE(6) 
replaces NTAPE(3) for EP-5, 6. If IOP = 2 in EP-5, 
several files are updated permanently; and NTAPE(6) 
should be saved and used as NTAPE(3) for subsequent 
runs. 

This is a scratch file required for EP-4. If IOP > 0 
for EP-5, it is required there as a scratch file also. 

This is an output tape containing only the system physical 
eigenvectors which are partitioned by substructure from 
EP-5. 

This output tape contains those selected system eigen- 
values and eigenvectors obtained for backsubstitution 
from EP-6. 

Special input tape for EP-1 when I OP = 2. Use previous 
NTAPE(2) data and new NTAPE(2) generated. NTAPE(7) 
may be used here if old NTAPE(2) copied onto this file. 

Special input tape for EP-2 when I OP = 2. Use previous 
NTAPE(4) data and new NTAPE(4) generated. NTAPE(7) 
may be used here if old NT APE (4) data copied onto this 
file. If NTAPE(7) used in EP-1 for NTAPE(IO), it may 
not be used again for EP-2 if the EP's are being executed 
consecutively. 

Special output tape for EP-5 when IOP = 2. This tape 
replaces NTAPE(2) for subsequent runs. 
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4. 2. 2 Data Locations on Tapes 


Location Matrix 

I 

II -- Basic data 


2 


1 

2 

3 

4 


3 1 


4 1 

2 
3 


5 1 

2 

3 

4 

5 

6 

7 

8 


ID DOF 
LAMDA 
PHYSUB 
IPDOF 


PHYPRT 


CPL1 A 
CPL22 V 
CPL12J 


K22 

K12-ll 
K 12- 2 ) 
K 1 1 

LAMSYS 
GAMK 
GAMR 1 1 

GAMR 1 2 


Vector for reordering substructure DOF's 

Substructure eigenvalues 

Substructure eigenvectors 

Vector for selecting substructure DOF's 

to be printed in physicals. 

Set repeated for as many substructures 
read into program. 

Partitioned set of substructure vectors 
assembled by IPDOF. 

Repeated in blocks of 100 DOF until finished. 
Partitioned coupling spring, CPL. 


Set repeated for each spring read into program. 

Partitioned stiffness matrix for reduced 
coordinates 

Partitioned stiffness matrix coupling terms 

Partitioned stiffness matrix for kept 

coordinates 

System eigenvalues 

Participation factors for kept coordinates 
Participation factors for reduced coordin- 
ates in partitioned form. 


6 


1 

2 

3 

4 

5 


M* 

2K12*Rll' 
2K12*R12 „ 
KZ2*R11 
K22*R12 j 


Reduced mass matrix 
Partitioned multiplications 
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Location 


Matrix 


7 

1 

R 1 1 

R matrix for dynamic transformation 
in partitioned form. 


2 

R 1 2 


8 

1 

PHYSYS 

System Eigenvectors 


• • • 


One record for each substructure. 

9 

1 

lambks 

Eigenvalues from backsubstitution 


2 

PHYBKS 

Eigenvectors from backsubstitution 


4 • • 


Repeat set for as many modes selected 
in groups of 30. 

10 


Optional input tape, NTAPE(2) 

11 


Optional input tape, NT A PE (4) 


12 


Optional output tape, NTAPE(2) 



FORTRAN LISTING OF DAMUS AND FLOW CHARTS 



This appendix contains the Fortran listing of DAMUS. The 
main program has been divided into four links with each link beginning 
with one of the six entry points. The main program calls the first link 
with all links being called in sequential order. 


Link 

Entry 

Subroutine 

No. 

Point 

Name 

1 

1 

ST Cl 

2 

2 

STC2 

3 

4 

STC4 

4 

6 

STC6 


The flow chart following the program listing shows the major 
subroutines called by the main program links. If the subroutine called 
requires input data, a flow chart for that particular subroutine has been 
included. 


A - 1 



D 1 0 6 = 06-73 _15.*5-4Z 


_STSL_MAJ?J. LINK 


CSTC STC HAI'J LINK 

COXM ON i F T R I E X ! T f f o U , NUT 1 , \ U T 2 , N'J T 3 , NUT 4 , N U f 5 7 M U T 6 , N U T 7 # U 1 J T 5 

1 , MUT9.MJT10; \"JT11, NUT12 

COMMON “/ML INEZ/ DUMYl /|_STAKT/ DUHY2<30) 

C Q M M 0 N _/ W 0 7 < V C / L N G 1 H W » W ( ^ q q ) 

""LNGTHW s 20 3 “ - 

CALL LjNK_ _ (_6HEMTRY1J 

S T 0? 

END.. .„ 


23 7 4 7 WORDS OF ME M 0 R Y_ US E D BYTHIS COMP IL AT 10 N 


A- 2 



5 11ELN g?S _C 0 UP LING __E N T_RY__1 


5 c 06 -13 15 ..54-2. 


CSTC 1 _ _ STIFFNESS COUPLING iENTRYJ, 

subroutine srci •- 

_J)IMEMSIOM_ A (103 , lOJ ) /SC 130,100 ) ,PHY(10 0, 100 ) 

1 , "ID (Yoc ) . I V ( r 00 > ; jvr t 00 ) T t C . 12 (Too > I T3 ( jWY 

_ 2 V _ _NNDOr (20,6)#^EEP(2»20)|KCPL(210J,KSK(20)*KSR(20) 

3 , I D P H Y ( 2 1 ) V I D C P L ( 1 0 0 ) 

COMMON! lENTR, IEX lT f FOD,MUTl,NUT2,NUT3,NUTj l N‘JT5 f NUT6*NUT7,M'JT.3 

I’, M‘JT9,N'jTlO l ..\iUTll,NUTl2 

COMMON / .* J L I u r: Z ^ QUMY 1 /LSTART/_DUMY2(_3q j 

" COMMON >.s’QPKVC/' LNOTHW, W(?OC > 

Data Kni,KD2/<D3,K34,KD5,KD6 f Kp7 l KD8/100.. 200;2*20,300, 21,6,100/ 

DATA L' 1»LL 2 'U 3;l> 4 *Ll 5>LL 6 »LL 7 , LL 9 »LLl°»LLll*LL 1 '^ 11 * 1,J ^ rj/ 

DAjA PHD<l>/i20HPHYl PH'12 PHY3 PHY4 PhY5 PHY 6 PHY 7 PRYs PH 
$Y9 ^HYl 0 PHY 11 PHY12 PHY13 PHY14 PHY15 PHY16 P H Y 1 7 PHY13 PHY19 PH 
5 Y 2 J / 

DATA ’ N SU 9 T 7 NO P L T , N NO 0 f\ I DP H V , I DC P L#"K C PL , K EE P V K 5 K , K S R / 5 3 3 « 0 / 

1JC0 FORMAT (1615) _ 

1 0 C 1 FOR M A T ( E i 0 Vo ) " ' 

2 0 C 0 FORM a T (///lOX'POD =',1PE12,4) _ 

2001 FORMAT ( / / / 1 3 X i 7 H I iij P U T T $PE NOS" ’=,1215 » /10 X17HUN I T A^S I Gf\|ED 

, 1, _ . 1 2 I 5 ) 

20 C 2 FORMAT </,//10X' IGNTR = 1 5 , 1 0 X ' I EX I T = 'I5> 

2003 PORMAT ( ///10X « ENTRY POINT' 1 5 , 4 X ' H A S BEEN COMPLETED.') 

2004 FORMAT < ///1QX ! I ENTRY' 5 ' I 5* 10X • NSUBS ='I5,10X'IOP =*I5)' 


C ; . 

1 CALL START 

READ (5,1001) FDD _ 

RE A iO ( 5 1 1 0 0 0 I N.jT 3 » NijTl , NyT6 » Ny T2 1 NjjT3 # Ny T 4 ,’NijT5 1 Ny T7 i Ny T9 
1, MUT10 i NUTll , NUT12 

READ '(5,1000) IrNTR.IEXIT 

WRITE ( $ » 2 0 0 0 ) FOD 

"WRITE (6, 2001 )“"(■ I , I = i", 12 ) , N JT3 , NUTl , NUT6 , NUT2fNUT3 , NUT4 , NUT5 , NUT7 

1 , :] : JT9,.N'JTio,wUTii, N UTi2 _ 

WRITE ( 6/2002 ) IEnTR, I5XIT 
IF (IENTR.GT.1) GO To 9002 
" ‘ IF ( IEMTR.E0.1) QJ TO 9 0 01 
REWIND NUT 3 

"READ ( N'JT 3 ) N S U R T , N N D 0 F , N CP L T , I D C P L , K C P L I DPHY , KEEP, KSK 

A , K S I? < M i M » M i * M 2 * M 1 P 1 » N S Y M » M S ¥ M 

c " 

9001 IENjR y s 1 

c 

C get sUB^TrUCTUrE MODAL DATA FrOM MASS AND sTlf/NEsS MATRICES 

c 

_ READ (5,1003)_ Ns:.jB s , IOP _ __ 

" CALL PAGE HD ' ~ 

WRITE (6.2004 ) I ENTRY » NSUBSj IOP 
IF ( I OP • EO • 2 ) GO TO 120 

CALL SU3DAT ( I QP . A , S , PHY , N.NDOF « I DPHY, I D , I V # JV , T ± . 7 2 • T 3 

~'1‘7' ~ NSUbS#N3U*3T,KDl , K!)4 , KOI , KD6 , LL* . NUTl , LL& . NJUT6 , FDD , PHD ) 

no T) 12 5 

"120 “NSUBT = 0 “ 




= 06-73 U5_, 54 7 STIFFNESS COUPLING ENTHy :.l 

CALL _SU3D_AT_ ( I OP A , S , PH Y , NNDOF , IDPHY, ID, !ViJV,T1#T2,T3 

1," "" ” N3UBS.;JSUBT < KiJi»K:j4,KDl,KU6,LL' 1 *NUT4,o,O.FOD, ai O) 

CALL OELSUO ( A , PHY , .'KJUOF , I DPHY, ! D , I V , JV , NSURS# NSUQT » KDi » KD? > O 7 , 4 
iV’" ’’ ~ L L i i NO T 1 , L L “ » N U T ° # L L 1 3 1 N U T 1 Q , L L 4 N U T 4 , N U T 8 ) 

125 REWI'JD SJIJT3 

~ W R I T c ' ( N U T 8 j'~ N S U R T » N N D 0 F , N C P L T , IDCPL , KCPLV IDPHy, KEEP , KSK ~ 

li K*P,n»M»m1*m2»m1P1iN^ y m»m^ y M 

CALL EOF 3 ( LLli N«JTl . LLn , NUT6 , LL3 , MUT3 > 

CAL L P P l NT I_ (_3i HOUR STRUCTURES READ I NTQ PRO GR AM* 6 , NNDQF (1,6) 

^ # NSljy T » 1 » 0 ) 

IF ( IEXIT.EQ. IFNTRY)*G0 TO 9990 


.9 g 0.2 _ CAL L J. rU ( 6 H E N T h Y 2 > 

9990 CALL PA OF HD 

W RITE < 6 ,2 0 03) I E N T R Y " 

STOP 

'END 


.^l^ORDS OF" HEilORV USED - ' BY~ THIS COMPILATION ' 



A-4 



6 = Q 6 -_7 3 15.-4-5.54. 


-SJ-C.Z. 


CSTC2 STC2 

SURROUT I ME STC2 


D I HERS ! ON A (_10Q # 10 0 ) ; S C10O , 10 0 ) , PH Y < 1 C 0 , 1 OOJ 

^ 1 2 ( 1 0 0 * 2 0 0 j i S 3 (2 0 1 0 0 ) 

I_V(100)»JV(100),KKPT(100»2).KRED(100»2) ,MRED<2) »MKPT<2) 

N.jDnr ( 20 » 6 ) * KEEP ( 2# 20 } « KCP|_ ( 210 ) , KSK C 20 j , KSR < 20 ) 

IDPHY( 21 > » IDCPL (100 > 

C0HM0TOENTR7IF.X [ T ,>'0D* NUT1 , NUT27 NUT3 , ‘N0T4VNCIT5 /NUT6 , NUT7 , NUTS 
NUT9 , NUT10 , NUTll, NUT12 
COMMON' /’YL INEZ/ DUMYi VLSTaRT/ QUH Y2 ( 3 0 )" 

COMMON / W 0 R < V C / L^R ThW , W < p n Q ) 

EQU I VALENCE ”< SS ( 1 ) , S < 171 ) > , < SS ( i >TS12 < I Vi )) . ( SS <1 0 001 ) . PH Y ( l", 1 ) ) 
DATA KD1, KD2. KD3, KIM .KD5,KD6,KD7,KD8/100720 0/2.20, 30 0,21, 6, 10 0/ 

DA T A LLl » LL2 » LL3 LLi # LL5 I LL6 , LL7 i LL9 i LLl 0 , LLH » LL12/11 *10 Q 0 / 
FORMAT (1615) 

FORMAT' (/ / /- OX * E N T R Y ’ PQ J NT * 1 5 *4 X ' HA'S BEEN COMPLETED, * ) 

FORMAT <///lOX' IFNTRY = » 1 5 , 10X ' NcPLS s » I 5, 10X • I_OP ='l5) 


10 00 
2 DO 3 
2006 

2007 VuR:' 1 A T '"(//7l0AC‘flENTRY ? i 1 5 > 

REWIND N WTO 

READ ( NUTS') ” N 70 fj'T 7 NN D D V , N C PL T , T D C P L,' K C P L,' I D P H Y , K E E P , K S K 

„ , KSR*ri»MiM- #M'j»M 1 P, » NSYM, MSYM 

l G0"'TO"<90D27o002. 9003^9004, $004 ,9004 >7 1ENTR~'~ 


C 

90 C 2 ' ~ ie;jtr y '= 2 

C 

C GE T I ND ] v I DUAL K* CONT W ItLUTruN s 'F;,OM COUPLING S P R ING S 
C 

1 READ ('571000') NCPl'sViOP' 

CALL PAGE HD 

W R I T E ( 6 , 2 0 0 6 ) » I E NT R Y ,N C P L S , IOP 
IF ( I0P.E0.2) GO TO 210 

CALL CPLSPG (A,S,PH v ,NND(1F.KCPl'VIDCPL,MCPlS;ncPLT,KD 1 ,KD 1 iKD4 
1, LL2,NUT2;lL1»NUT1) 

I F ( I E X I T . E 07 I F N T H Y . UR 7 "j 0 P , fc Q V ±T ~G 0 ' T 0 ? 2 0 

GO TO 23 0 1 

21 0 " McPLT = 0 

CALL CPLSPG ( A , S , PH Y , NNDOF , KCPL , I DCPL, NCPL.S [>UCPLT r KD.1 , KD1 , KD4 

' i7" “ LL4 1 NUT/J , LLi , NUTi ) 

CA LL DE 1 CP 1 (A,Nf|DoI , IDCPt ,KCPi , lV |N CP L S ( ‘jCPi T,KDl,KD2,KD7,KD8,KD 

u LL 2 » NDT 2 , LLll i NUTil j LLA » NUT 4 , NUTll J 

220 REWIND NUT 3 

; WRI f E ( N U T 8’7 N S 0 8 T 7N N D 0 F ,~M C P L T , IDCPL7 KCPl.I OPHY/KEEP, KSK 

1 KSR»N»M»Mi»-H2*MiP 1 iNSVM«MSYM 

— CALT"E0F3 ai'8,NUT870-,0;070) 


?30 CONTINUE 

- Ca LL E0F 3 “( LL2 ; HUT2 , 07070 , 0 ) ‘ 

C Ait CNNCiV t A.NnDOF,KCP, , KD4 , KD1 , KD4) 
IF ( I EX I T. EQ . TENTH Y ) GO 'TO 9990 ‘ 


C 

9 00 3 rENTR Y “="3 

C 

C K* AssEMBl.Y 


A-5 



OOO OOO 0.0 0 


>» 06-73. 15*554 


..SIC 2 


c _ 

CALL PAGKHD ~ ' r " " ' ' “ ~ 

WRITE: ( 6 » 2 0 C 7 ) IFNTRy 

CALL READ 1 M ( KFF.P • NH , MXSU8 . KD3 . ) ‘ 

CALL ! NDCFS ( KEEP* MXSU3* KSK, KSR, N, M, NSMj. , Mi » M 2 * MiPi» NSYM» MSYM,_KD3 ) 

REWIND NUT* ~ " " 

WRITE (NUTS) \'SURT f NND:)F,NnPLT, IDCPL,KCPl. I DPH Y , KEEP , KSK 
1 ,' ' " K 3 R * *1 # M * fl 1 » ft 2 » M P x » N S Y M , M S Y M 

CAi.L-.E 0 r 3 (_ll3, nut 8, 0, 6,o,o_) 

: assemble k^p* 

CALL AsHRLl (2,SS»A,K h ED,M r ED,KEFP # NND0F,K S r;KCPL, I V , JV, M , MX s uB _ 

1, HSni,:iSYM,KUl,K'j3,KD4,i f UL3,NUT3’;LL2,NUT2,ULl,NUTl) 


Assemble kir* 

CALL' ASMBL2 < 312 , A , KKPT , MKPT , KRED , MRED , KEEP , K S'kT'K S R K C PL ,' j V , J V NV M 
1 V _ M 1 . M 2 , r M 1 P 1 . M * SUB , ;-j S M 1 , k D 1 , K D 3 , 2 , L L 3 , N 1 U T 3, LL 2 f N | U T 2 )_ 

assemble .Kii* 

CALL A S H B L 1 (l.SS, A , KKPT , MKPT , KEEP, NNDOF , KSK'.'KCPL , I V , JV , N , MX S UB 
1, NS l< l,.NSYM,KDl,KD3,Kn4 < 4,LL3,NUT3VLL2,NUT2 # LLl l NUTl ) 

CALL E OF 3 ( LL3 , NUT 3 , 0 , 0 , 0 , 0 ) 

IF (JKXIT.eq. IEMTRY) GO TJ 999 q 

C _ 

■ 9yU4 CALL LINK ( 6HEM T RY4) 

C .. . „ 

9990 CALL P AGE HI) 

WRITE (6,2033) IFNTRY __ 

STOP ' ' ' ’ ~ 

END 


WORDS OF MEMORY USEn BY THIS COMPILATION 


1 


A-6 



•Q 6-7.3 15^55.a 


ST1F FNESS__C0UPL I N.G EN.T.BX-4. 


CSJC4 _ ST IFFNESS .COUPLING ENTRY 4 _ 

‘SUBROUTINE ST C 4 ‘ ” 

DIMENSION A(1C0,100) # SC100,100) , PHY (100, 10 0 >;si2 (10 0,200} 
”S?:iT200. 1 QOi»R(200.‘>0).R2t2C0.50).SSC 2 0 1 00 > 

2 V I'0(100>_. IV(100),JV(10 0).T1(100)»T2(103>#T3<100> 

3 , NMnoF ( 20 •. 6 ) . KEEP < 2 ,‘20 ) , KCP U ( 210 ) ; KSk C 20 > I KSR { 20 > 

4, LJPHY(P1),IDCPL(100),PHD(20)_ 

5, LaNSYS(IOO) 

COMNOJi JF.NTR, l E X I t *KOD»NUTi »NUT 2 ,NUT3,NUT4.NUT5.NUT6»_NUT7, NUT8 

1, NUT9*NUTlO,NUTll ( NUTl2 

DO MM ON /ML INEZ/ D U M Y 1 YLSTaRT/ DUMY2(30) 

CO M M ON/ W 0 R K Y C / L N G T H W » H ( ^ g g ) 

FQU IvALcMCE <SS(1> *S( 1»1) ) f (SS(1) #S12( 1|1J) » <SS{1) #S2ltl|l) ) 

1, " < A < 1 » 1 > # 1 1 1 » 1 ) ) 

2, iSSjlOOO.1 J , pH_Y(i,i) ) , (SS<.100.01).. t R2£l i .lL ) 

real” LAMS yS , LANDAU 

DATA KD1 < KD2.KD3,KD4,KD5 i KD6,KD7,KD8/100, 200V2, 20 , 30 0, 21 , 6 , 1 00/ 

DATA" LL1.LL2, LL3.LL4 , LIT , LL6 , LL7 , LLP * LLl 0 , l.L U , LLl2/U*l 0 0 0 / 

DATA PH1K1 ) /120HPHY1 PHY2 PH Y3 PHY4 P H Y !> PHY6 PHY7 PHY8 PH 

Sy 9 PHY 10 PHY PHY l2 ' PHYi3 "PHYi4 PHY'i5"pHYi6 PHYi? PHYi® PHYi^PH 

$Y20 /_ . .. 

lOCO FORMAT (16155 

2 nr 3 rORMAT ( ///in>: ’ENTRY PO HIT • I 5 , 4 X ’ HA S BEEN COMPLETED,') 

2oC5 FORMAT ( ///10Y ’ JENTRY ='15,10X'MD1 s ' 15, 10X »MD2 =*I5*10X»IOP = 'I5) 

20C7 FORf-iAT ( ///IGX ' lEfjTHY = 'I5) _ 

' REWIND iJ'JT 8 

READ (NUTB) NSUuT;nnDOK,NCPLT, IDCPLiKCPL, idphy,keep,ksk__ 

, , KSR » N » M > M i » M'2 » M j.P - » N5 YM » MS YM 

no TO (9004, 9304;9004 f 9004, 9005,9006), I ENTR 

c " 

90 C 4 IF^T rY __ s 4 — _ 

C 

c fo h m r fo r _dynamic transformation 

c 

CALL PAGEHD __ _ 

WRITE (6*2007) I f : N T R Y 

CALL READ (LAMDAOj NH. NO, 1» 1 ) _ 

CALL R T A P S S ( S S » M R * NC * A N M * i , L L 3 * N U T 3 ) ~ ' “K; 

CALL RM A T _ (SS|R*LAMUA0*M.KD 2 , 1 *LL5«NUT5*2»LL3#NUT3)_ _ 

c * 

CM. 

C 

CALL MSTAR U,s21*!'l,M,MlPl,KDl,KD2,l,LL4,MlJT4,l.LL5,NUT5) 

C 

C K* _ . 

c 

CALL KcTAo (sl27p*s<:f s» A, p 2»N,M,M1P1,K01;KD?,2»LL4,Nut4 

I ,' 2 . LL3 , NUT 3.1, LIT , NUT 5 , 1 ,LL3, NUT3) 

C _ 

C “ s'OLVE' F Or EIGENVALUES AND EIGENVECTORS 

c 

CALL' E IGVEC "( A ; "s ; S 2 1 L A M s Y S . T 2 , T 3 , N , M , M 1 , M 2 , M 1 F> 1 , K D 1 , K D 2 



6 = 0 £> r.7. 3 15^55.3 S 1 1 F f_N E S S . .C O.y P L I f I G JEN.T « Y A 

1 . _ 5 * L L 3 , f4 U T 3 , 1 , L L 4 ( N U T 4 , 1 , i. L 5 N u T 5 F 0 D > 

CALL E0F3 (LL3#N'JT3»ri#Q#g*n) 

IF ( IEXIT.EQ. IENTRY) GO TO 999 n 

'c * “ ' ' ‘ 

9 J 05 IENtRY = 5 

: c ' ~ ----- ■ — - “ 

rEAD <5.1000) M D 1 » M D 2 * I 0 P 

CALL PAGEHO ' ~ 

WRITE (6,2005) IENTRY, MD1, ND2, JOP 
I F ( I OP . E Q . 0 ) GO T 0 52o 

CALL DELPHY (A,PflY,NNDOF» IDPHY,T t , ID# IV, jV»N5UnS#KD 1 .KD4iKD6 

1,' 'LL4 ( N'JT4,UL5,NUT5,LL1#NUT1) 

IF ( IOP.FQ. i ) GO TO 5,0 

C A LL ■ UPDPHY *r A . NND'jr i 3V, fjsu a s # NSU B T # KDi #'K‘n4 ,'LL12TNUT12 , LLl #"NUTi“ 
1. I l 5 , M u ~5) 

' RF.wTfTn 1 JUT 8 “• ' 

WRITE (MIJT3) 'ISURT.NNDOF.NCPLT, IUCPL.KCPL, I DPH y , KEEP , KSK 

1, " K S R , M # M , M 1 , M s , M 1 P 1 , N S Y M , M S Y M 

CALL EOF3 (LL4 , NIJT4, LL'3# NUT3, LLl«-, NLJT12 ) 

“"5,0 LLo = LI. 4 “ 

A NUT6 = MUT4 

c " 

C GET PH y 5 ICA!. S 

c 

520 CON T !i\'uE 

C ALL PH YSCL ( PH Y , $7 A , R i NNUOP , KFE P , KSK # KSR 7 I DPH Y > f\l'» MDl , MD2 , KDl 
1 1 __ K02#KD4 ,kD 3,KU6, l,LL7#.NUT7,6. LL3 l NUT_3,tfL6,_NUT6,PHD>_ 

C aLL "EOr3 " ■■■(LL7,.NUT7, 0 , Q , 0 # D ) 

IF < IEXjT.F/i. IgNTRY) qO TO 9990 

: c 

9006 CALL LlNK (6 hE’JTRy6) 

,’c ‘ " ’ 

9990 CALL PAGE HD 

'WRITE "< 6,2093 7 "IENTRY 

STOP ___ 

END" ' 


WORDS ’ OP HEifCR Y USE j' B Y“'TH I S ' COMP I L aT I ON 


6 • Q 6-73 i5. t S 62 SJC 6_ 


CSTC* _ STC6 _ _ 

SUBROUTINE STC6 

D I M E * 1 3 I 0 f J A £ 100,10.)) .SC100.100 ) f PHY (100 ,10 0 )VS12C 100,200) 

""s'>ir20n/ 1 PO)»GAMA( 300»30)»SS(2b 1 nO) ~ ' 

2 ; IO( 100 ) f Tv (100 ) . Jv ( 100 ) ,T1 (100 ) ,T2(10C) ,T3( 100 ) 

3, ' NNno( r (PC,6),KEEP(2.20) ,KCP L <2lO),KSK(20).KSR(2n) 

4, ' IDPwv(2l 5 , IDCPU(IOO) , PHD (20 5 

5 ; " •LAMSYS(l'.lO)- ' 

COMMON 1ENT1, l EX I T , FQD > NUT < , NUT? # NUT 3 # NUT 4 #' MUT5 # NUT 6 # NUT7 , NUT 8 

'*T, ~ “ " M Ll T 9 y.NUTlO , NUT11, NUT12 ” ~ 

• COMMON /N LIMP:// PUflYl /LSTART/ DUMY2(30) 

• COMMON' /MOR-;'.'C/' L NGT H Vii'W( 2 no) ‘ “ 

F 3 U I VALENCE (SS(l)r'S£l,n),£SSCl)|Sl2(l,in,(SS(l>,S2l(lU>) 

1 ,' ( A ( 1 . 1 ) I OA m A { 1 , 1 ) } 

2 ,' (SSllOOOii ,PHyJ 1 i 1 ) ) | 

REAL LAM SYS 

DATA KD1# KD2. KD3; KD«; KD5, KD 6 , KD7.KD8 /100;20 0 y 2,20, 300, 21, ft, 10 0/ 

DATA''LL1.LL2, LL3,LL< , LL5 , LL 6 , LL7 , LL9 , LL1 0 , LL11 . LLl2/ll*10L0/ 

data PHD(1 )/120HPHYl PHY2 PHY3 PHY< PHY5 PHY 6 PH Y7 PHY 8 PH 

i $ Y 9 P H Y -1 0 PH YL^ PH Y 1 2 PHY 1 3"PHY 1 4 PHY^5' PHY '^6 PH Y^7 PHYj.0 PHYj.9 ph 

*Y20 / _ X __ ___ _ _ 

1 10 00 F f j RM A T (loi?')” 

2003 FORMAT ( /// 10 X 'ENTRY PO I NT • I 5 , 4 X ' H A S BEEN COMPLETED,') 

2UC8 FORMAT ( / / / 1 0 X ' IENTHY =' l5, 10x ' NBKSB =' 1 5 j 

REWIND NIJT9 _ 

Li. 9 V ' 1 ~ 

reWjNd nuts 

READ (Nut 8) ’'•) S ij 3 7 , N N D 0 F , N C P L j I D C P L "# K C P L , ID P H Y ", K E E P , K S K 

1 , KSR , N , M , Ml , M2 , Ml PI , NS YM , MS YM 

c " ' " 

9-J06 I ENtRY = 6 

C ' 

C BACKsUDsTlTUlION _ 

C 

READ (5,1000 ) N o K s B _ _ _ *1 

CALL PA SEMI) " " ” ' 

WRITE (6,2003) I ENTRY, NBKSB 

CALL RTAPDS ' (LAM5YS,NR.NG;aNMV1|5,LL^'.NUT3) 

DO 49 Q J .) = -! » NBKSB 

CALL' BaKSUB (G. 4 MA,S: 2 ,S 2 l,S 3 ;S,LAHSYS,lD,rV,'JV,Tl,T 2 ,T 3 ,N,M,NMn 

1 , KT5,K i )1,KD2,6»li 3,mHT3 i 2,ll3.N ut 3,1,lL 3 » nUT;5) 

CALL WTAPPS' (l 2 il»MMn, 6 WLAMMKS,i ,LL9,LL9,MUT9) 

CALL PHY3CL (°HY,S . A , 3AMA , N.'JDHF , KEEP ,-KSK , KSR ,' I DPHY » N , l.NMl) 

IV" KD17KD'5,KD4,KD 3 iKD6.LL9,LL9,NUT9;o,0,O.LL6 i ;hiT6,PHD) 

49 0 CONTINUE 

■ ' ca 1 - 1 - E 0 F 3 . { LL 9 , not q ; 0 1 o7o » o ) 

9990 CAi [ PA-3EHD 

y R I T E (6,2003 ) I E N T R y " " ” 

STOP 

END ' ~ ‘ 

WORDS OF MEMORY 'USED BY THIS COMPILATION " 
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RF = Read formatted input data. 


RM = Read matrix data using 
FORMA routines. 
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CPLSPG 


ISUB1, ISUB2 
NS1, NS2 


Coupling spring 
stiffness 


Assemble K ^ 

Assemble 

Assemble 


Return 


W T 

L k cpJ L 0 

W T 

OcpJtAl 

W T 

[ K CP J 1*2 1 









RM 


RM 


Vector containing 
substructure ID's 
for print changes 


Vector for 
selecting DOF's 
(one per sub. ) 


RM 


Vector to 
select modes 




Back subs titutio: 














APPENDIX B 


SUBROUTINE EXPLANATIONS 



This appendix contains a brief description of the subroutines 
specifically generated for DAMUS. Some of the subroutines may be 
used in a general FORMA program. The following list of general sub- 
routines were used by DAMUS. Detailed explanations of each routine 
are contained in the subroutine comment cards. 


Subroutine Names 


ALPHAA 
BT ABA 
BTABZ 
BTAZ 
BTB 

COLMLT 

COPY 

DCOM1 

DIAG 

EIGN1 

INTAPE 

INV4 


LT APE 

MULTZ 

NEGATS 

NEWMOD 

ORDER 

PAGEHD 

PRINT 

PR INTI 

READ 

READIM 

REYSYM 

RTAPE 


START 

SYMLH 

TRANS 

WRITE 

WRITIM 

WTAPE 

ZERO 

ZZBOMB 


B-l 



Subroutines for DAMUS 


ADDLAM 

ASMBL1 

ASMBL2 

ASSEM1 

ASSEM2 

BAKSUB 

BSOLVS 

BTABSA 

CNNCTV 

CPLSPG 

DCMSYM 

DE JLCPL 
DE LPHY 

DELSUB 
EIGYEC 
EOF 3 
GET PHY 


Adds substructure eigenvalues to diagonal elements 
of single-subscripted, symmetrically stored matrix. 


Assembles K^, K partitions in symmetric form. 

Assembles partition, double subscripts 

Adds K contributions into correct locations of K , K 

Jr 11 lZ 

Adds K„ contributions into correct locations of K 

-L* 12 

Backsubstitute s eigenvectors 

Solution for Z of equation [a]] £ Z~\ - [b^ where 

[Aj is symmetrically stored. 

T 


Matrix triple product Z = B A B where A is symmetrically 
stored. Z is double- subscripted but only upper-half is 
returned. 


Calculate connectivity for substructures. 

Generates spring contributions of 2 substructures in 
partitioned form. 

Decomposes symmetrically stored matrix A into factors 
where A = L D L/^ 


Changes data files when changing a coupling spring 

Changes data files when changing physical vectors to 
be printed 

Changes data files when changing basic substructure data 

Gets eigenvalues and complete set of participation factors 

Write end of file for up to 3 files at a time 

Used to generate substructure vector matrix for obtaining 
physicals 
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IDFILL 
INDCES 
KEEPV 
KSTAR 
MODE 2 
MSTAR 
MULTBS 
NEW PHY 

PHYSCL 

REVA2Z 

RMAT 

ROTATA 

RTAPDS 

RTAPSS 

SUBDAT 

TRANSA 

UNPAKS 

UPDFIL 
U PD PHY 
WTAPDS 
WTAPSS 


Special routine which fills up vector for omitted values 
Generates initial data and index locations 
Generates vectors for assembling K matrix 
Generates reduced stiffness matrix 
Eigenvalue /vector routine with mass matrix options 
Generates reduced mass matrix 

Multiplies BZ = A * BZ where A is symmetrically stored 

Generates new set of eigenvectors in partitioned form 
for getting physicals 

Calculates physical eigenvectors 

Revises matrix A into matrix Z, where Z may be a single 
or double subscripted matrix 

Calculates for dynamic transformation 

T 

Special triple product Z = R A R when a 3X3 submatrix 
of direction cosines is input for R. 

Special matrix read tape routine - double subscript 

Special matrix read tape routine - single subscript 

Generates substructure eigenvalues and vectors from 
basic substructure mass and stiffness matrices 

T 

B = A where A, B occupy the same core location 

Unpacks symmetric single- subscripted array into 
symmetric double subscripted storage (upper half only) 

Updates data from two files to one in a predescribed order 

Updates data file for physical vectors 

Special matrix write tape routine - double subscripts 

Special matrix write tape routine - single subscripts 
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APPENDIX C 


SAMPLE PROBLEM 



This appendix contains a sample problem using all 6 entry points 
of DAMUS. 

The sample problem used to illustrate the use of DAMUS was 
Problem 1, a 20 DOF longitudinal rod model consisting of 2 substructures 
(10 DOF/ substructure). The output data printed by DAMUS follows a 
listing of the actual data cards used to execute the program. 

Since this was only a test case, no output tapes were saved and 
the same logical tape unit was assigned to the input files 7, 8 and 9. 
Normal execution of DAMUS would require different logical unit numbers. 


C- 1 



INPUT DATA CARDS 


SICA 0 1 E..I. MJHAK 

stiffness coupling CHECK case FROM ENTRY POINT 1 
2 U nni LONGITUDINAL. ROD 2 SUBSTRUCTURES 


u . 


1 0 

1 1 

12 13 

l 

6 


2 



1 

1 0 

2 

I D D OF 

1 

1 

1 

1 

-1 

noon n o o u o u 


K1 

. 10 

10 

1 

1 

200 0 . 

2 

1 

-2000. 

3 

2 

- 2 ono . 

4 

3 

-20 011. 

S 

4 

-2000. 

6 


-2000 . 

7 

6 

-2000. 

8 

7 

-2000 . 

9 

8. 

-2000. 

10 

9 

-2000. 

(I 0 0 0 0 0 0 0 0 o 


Ml 

1 

10 

1 

1 

1 . 

1 

5 

1 . 

1 

9 

1 . 

oo oooun ono 


i p nnF 

1 

3 

i 

1 

1 0 

n o n 0 n i) n o o o 


9 

c 

1 0 

2 

I U D U F 

1 

1 

1 

1 

-1 

n o o n o o n n o n 


K 2 

1 0 

1 0 

1 

1 

20 0 0 . 

2 

1 

-2000. 

3 

2 

-2000. 

4 

3 

-2 0 0 0 . 


1 4 1 !> 16 16 


OROERS DOF’S FOR 


STIFFNLSS MATRIX 
-2000 . ' 

4 0 0 0 . 

4 0 0 0 . 

4 0 0 0 . 

4 0 0 0 . 

4 0 0 0 , 

4 0 0 0. 

4 0 0 0 . 

4 0 0 0 . 

6 0 0 0 . 

MASS MATRIX 

1 . 

1 . 

1 . 

SELECTS DOF TO BE 

-10 


ORDERS DOF’S FOR 


STIFFNESS MATRI X 
-20 0 0 . 

4 0 0 0 . 

4 0 0 0 . 

4 0 0 (I . 


16 - 


EP- 

ATTACHMGNT points 


-2000. 
-2 0 0 0 . 
-2000. 
-2000. 
-2000. 
-2000. 
-2 0 0 0 . 
-20 0 0 . 


1 . 
1. 


PR I NTE 0 


attachment POINTS 


-20 0 0 . 
-20 00 . 
-20 0 0 . 


N. B. * In Column 71 cancels call to Forma subroutine PAGEHD. 
• * C-2 



5 

A 

-2000. 


4 0 00 . 

-2000. 


6 

5 

-2000. 


4 o n n . 

-2000 . 


7 

6 

- 2 0 0 U . 


40 0 0 . 

-2000. 


8 

7 

-2000. ' 


4 n n o . 

-2000 . 


9 

8 

-2000. 


4oon . 

-2000. 


1 0 

9 

-2000. 


2 0 0 0 . 



0 0 0 IJ n () 0 1) 0 0 






M2 

1 

10 


MASS . MA TR I X 

- 


l 

1 

1 . 


1 . 

1 . 

1. 

l 

5 

1 . 


■1 • 

1 ♦ 

1. 

1 

9 

1 . 


1. 



non n noon no 






I poor 

1 

3 


SELECTS DOF TO 

BE PRINTED 


i 

1 

1 n 

-10 




ojj o u o u o u n n 






l 

l 

1 

1 3 

1 

1 


EP-Z 

KC PL 

2 

2 


COUPLING SPRING 

ST I FFNESS 

MATRIX 

1 

1 

2 o n n . 


-2nrn . 



2 

1 

-2000. 


2 0 0 0 . 


’ 

n n n ti n n n ti n o 






KEEP 

2 

2 


SELECTS KEPT/REDUCED DOF'S 

BP-l 

t 

1 

2 2 




2 

1 

8 8 





non n o o n o n o 






L A M 1) A 0 

1 

1 


P* *2 REDUCTION 

FREQUENCY 


1 

1 

0 . 




n o n n n 0 n u o n 






i 

A 





EP-S 

i 

I mode 

1 

4 


SELECTS MODES FOR BACKSUBST'] TUTI ON E?-b 

i 

1 

1 U 

- A 




ouunououoo 













j 

o 


Ui 

UJ 

CL 


O 

CL 


113 
x 
x 
. u 


o 

CL 


cr 

o 


o o 

O O 


. I 


t CT 
•< 

X 
z> 
rose: 

H • 

K>-J 
U> - 
O LLI 

UI >- 
*— m 
«< 

Q 2 : 

3 

: a: 

1 X 


i n 


2 

oc 


o o 
o o 
o o 




cz 


X 





O O 











O O 




in . 


cr 

i 

• 1 . • 



O O 




. cr 


»— 

- 




• • ’ 

‘ 



LLI 


< 

a 

0 0 O a 0 0 0 



H H * 



c\ 

C 


2 

> 



X 





X 






•— » 



1 


• O 


in 




cr 







cn 


wntow«rorofo 


1 — 

O O j 


■ 




LU 


00000000 


<1 

0 0 ! 



ir 

1 i 


-z— 




2 

t 




J 


u. 


LL' LL) UJ LL 1 UJ LU LLI UJ 



LU UJ , 



»— 



u. 


OOOOOOOO 


cn 

O O 



CL 



*— 


COGOOOOO 


cn 

O O 


»— 

a 



V— 


doc OOOOO 


< 

O O , 


2 

2 

1 

, 

CO 


OOOOOOOO 


2 

O O 

, 



1 



- 

0 0 c 0 c 0 0 0 



O O 

1 

0 in 



• , 



OOOOOOOO 



a 0 


CL UJ 






OOOOOOOO 



0 0 


cr 


; 




cccccccc 



c c 


>- D 

« 




• 




• • • 


cr i — 


! 



0 

C\IC\IC\jC\J(\jC\)C\J(\i 





h- 0 


i ^ 

• 



K 1 1 1 1 1 1 • 





■ 2 X 


1 ! 

J 








* uj cr 


: *-< : 

• : 







, 



( 

• 


m 


1 



* 

' X CO 

a 

; 1 



CD 

oooooaooo 



OOC 3 


, 0 cr 

r-i 

1 1 









cr x 

. 

• X ; 



UJ 

LU LLi UJ UJ LL) UJ LU UJ UJ 



UJ UJ Ui 


Li. in 


< , 



0 

ocoooococr. 



croo 



I 




0 

O OOOOOOOO 


«'-*■ 

coo 


t LU CM f 11 




0 

OOOOOOOOO 



coo 


(/> 

| 



0 

0 

OOOOOOOOO 


0 

000 


. 

c 

, ; 


H 

0 

OOOOOOOOO 


rl 

coo 

• 

; u 

c 


’ 


0 

COCOC.OOOC 



000 

1 


2 

j 1 

1 1 


0 

cr ooocoooo 



000 

1 

: r 

1 a 

J ; 

' * 

X 

a 

OOOOOOOOO 


X 

oca 

' 

* u a 

ui i ! 

» 


• 




• • •- • 


; ui 0 


Ll. ; 


Q 

c\ 



rl 

1 

1 

1 

i 

' x cr 


O 1 

: ! 

H 

1 




! : 


■ u 


Q } 

1 j 

' 



i 


J 1 

1 

- 

J 

! O 





; 


: I 


C 5 <£ 

1 


1 


ro 

fofofowtofofoion 




i 

i ~5 

' 

I - 

1 

1 X 1 

1 

t , 



oaoooaoac 

i 


i 

q 0 q 

• 

; -i c 

| 




LU 



• 

lii LL* U_*' 


. CL Z 2 


i tr j 

1 1 


a 

OOOOOOOOO 



000 


1 => •- 

j 




0 

OOOOOOOOO 


1 

oca 


1 0 — • 

I 

: < ! 

1 


-a 

OOOOOOOOO 



000 


! O 13 

j 

i 2 Z j 

1 


0 

oaoooocoo 


H 

a 0 oi 


2 ? 



! 1 

iC 

a 

OOOOOOOOO 

• 

2 

0 0 cq 


i in a 


• cr I 



0 

c. coGooooe 


; 

00 a 


1 to _ 


i ui ; 


X 

0 

OOOOOOOOO 


X 

0 0 OI 

i 

: ui 

I 




0 

OOOOOOOOO 


•-» 

oca 


• 2 u_ 

i n 

! Ui I T-« 


cr 

- 



cr 


• 



1 »- 1 » 

' . • 

y~ 

C\J 

C\j 00 cv CM CV CN. CVJ CM CV 


i - ’ 

H r 4 t-i 

1 

■ U- a 

ca . z : 

2 



l 1 . I II It It- (fit II 


•<. 

j 1 

i 


1 X 



2 


: ; ! 

• • 

2 ; 

J 1 

* • 

!*- a 
j v. c\j 

. <n : 1 

— t • *- \ 

' Q 
• < 

*-* 


H (V fO ON CD’ O' 

Q 

►H 

«4 IA 0 >| 

a 

<c 


1 


Q 

cr 


LU 

cr 


: o 

2 


1 : i 

HCMr)'Tin'0Nt00'O 


C3 

cr 

< 

o 


: C-6 


. in 

cr 


: O 
: Q 


Q-! j 

S i 


Q 

cr 


a 

a 

2 

Ui 


H 


! I 





o 


j 


I 


I 


! 


I 


U1 

00 

o 

o tn 

*4 O 
r- 


Ui 

fN. 

— 'O 

O «H 

O 

^ ' • 

CO 


Ui 

o 


CD 

< 

a 


o o 
! in i 


Ov 

'O 


Ui 

C\J 

, O' CO 
! 1 TO 


i CC 
\ •< 

, X 
. ^ 

ro* 

r^ 

T-i • 

ro -> 

in • 

O LU 


T- CD 

. 

e: 


<< 

u 

►4 

1/5 

i 

£ 

2 

I 

2 

cr 


: 

; 

fO 


» .,H 


w4 

i 

i 

i O 


o 

; 

• O 

1 


i Ui 

! 

! uj 

1 

Ui 

i 


, V 


*r. 

i 

CM 



CO 

~ ro 


in 

j 

; co 

CM 

CO O 

1 co 

rl 


• 

i CO 

i 

CM 

! 

ro 

• 

• 

• » • 

w • « 

• ^ 

- ■ 

i 

1 

! 'O 


00 

* 1 

i 


i 

j 

i w 

i 

i ! 

1 *•*< 

1 

I ^ 


» 

o 

j 

J o 

i 

! o 

j 

. 

: UJ 

; 

• i 

UJ 

; 

’ Ui 

1 

■ 

c 

i 

! CM 


CC 



! -O 

— <3 

i 

ro 

i 


ri 

n- cm 


H 



CC 


o 


i CM 

1 

w 

' • 

w • 

>**■ 

• 


1 

> in 


rs 


I H 



TO 



\ 
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/ Headquarters: Valley Forge, Pennsylvania □ Daytona Beach, Fla. □ Cape Kennedy, Fla. 
□ Evendale, Ohio □ Huntsville, Ala. □ Bay St. Louis, Miss. □ Houston, Texas 
□ Sunnyvale, Calif. □ Roslyn;Va. □ Beltsville, Md. 
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